Level curvatures, spectral statistics and scaling for interacting particles. 
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The mobility of two interacting particles in a random potential is studied, using the sensitivity of 
their levels to a change of boundary conditions. The derealization in Hilbert space induced by the 
interaction of the two particle Fock states is shown to decrease the mobility in metals and to increase 
it in insulators. In contrast to the single particle case, the spectral rigidity is not directly related to 
the level curvature. Therefore, another curvature of topological origin is introduced, which defines 
the energy scale below which the spectrum has the universal Wigner-Dyson rigidity. 



PACS numbers: 72.15, 73.20 

The interplay between disorder and interactions is a 
central rapblem in quantum transport. Recently, it was 
proposed!!! that interactions could favor derealization in 
disordered insulators of large (one particle) localization 
length L\. This was first understood for the simple case 
of two interacting particles (TIP) in a random potential, 
where a fraction of the TIP-states do haveS a localization 
length L2 larger than L\. This result is surprising since in 
disordered metals, repulsive interactions are expected to 
reduce the conductivity^! and in disordered insulators, the 
singularity of the density of states which leads eventually 
to a gap, also reduces the quantum transport. This de- 
localization effect can be understood by considering the 
problem in the Fock space, where the many-body states 
are delocalized by the interaction. This is a very general 
idea that we use here in the restricted framework of the 
two-body problem. By Fock states, we simply mean the 
2x2 Slater determinants built out of the one particle 
states, assuming spinless fermions. These are eigenstates 
in the absence of interaction, and form a basis of the TIP- 
Hilbert space. This terminology of Fock states might not 
be the correct one for this two body problem, but we use 
it both for brievety and because it will become appropri- 
ate for the general N-body problem We consider bare 
particles, but the extension to dressed quasi-particles cre- 
ated from the Fermi vacuum is straightforward □. In the 
presence of interactions, characterized by the parameter 
U, the Fock states are broadenedtffl, and the TIP eigen- 
states have a finite projection over (typically) gi = r/A| 
nearby in energy Fock states. Here, T is the broadening, 
characterizing the local density of states in the Fock ba- 
sis (assuming a Breit Wigner form) and A2 the spacing 
between the Fock states effectively coupled by U. This 
derealization in the Fock basis characterizes both metals 
and insulators. But the derealization in the real space 
occurs only for insulators, since the Fock states are them- 
selves localized in real space. 

In order to understand this derealization, a scaling 
approach of the TIP-problem was proposed □ as a gen- 
eralization of the Thouless block scaling picture. For 
system sizes L\ < L < L2, a certain conductance 
g%{L) — T/Aj was defined and assumed to obey the 
usual Ohm's law (i.e. the one particle scaling theory 
of localization) in order to recover the derealization by 
interactions. In that work, as in an effective a model for- 
mulation developed^ later, L\ was implicitely assumed to 
be the smallest resolved scale. Therefore, the meaning of 



172 is unclear for L <C L± i.e. in. the metallic regime. 

Independently, an analysisQ of the statistical prop- 
erties of the TIP-spectrum in the metallic regime has 
shown that the universal Wigner Dyson rigidity occurs 
for g 2 = -%e consecutive levels, provided g% > 1. In the 
absence of interaction, the spectrum is essentially of the 
Poisson type for energies A2 < E < Ai where Ai is 
the one particle mean level spacing. For the one particle 
spectrum, one knows that Wigner-Dyson rigidity occurs 
for g\ consecutive levels, where g\ is the dimensionless 
conductance of the non interacting electron gas. For the 
TIP-case, gi plays the role of g\, as far as the spectral 
fluctuations are concerned. We will show in this letter 
that, though the TIP-spectrum becomes more rigid for 
stronger interactions, the transport in real space is re- 
duced when L < L\. 

To that purpose, we describe the transport in terms 
of the sensitivity of the levels to a change of boundary 
conditions a. This is equivalent to closing the system as 
a ring of length L pierced by a dimensionless Aharonovh 
-Bohm magnetic flux <f> = $/$o where $0 = hc/e. We 
define the two-particle level curvature C2 (E) at energy E 
as the zero flux curvature of the integrated two-particle 
density of states N(E, cf>): 
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An equivalent expression was shown □ to describe the 
d.c residual dimensionless conductance g\ of the non in- 
teracting electron gas under very general conditions set 
by the random matrix theory. It is therefore natural to 
try to extend this description to the case of interacting 
particles. 

The main results of this letter are as follows. AssuHiing 
a Breit- Wigner form for the local density of states M, we 
obtain for the two particle curvature C2 the expression 



C 2 (E) = C^(E)-f-I(^f 



(2) 



where C<^ '{E) is the value of the curvature in the ab- 
sence of interactions, B the kinetic energy scale (band- 
width) and /(f,f) w Sl p- for g x > 1 and of order 
— Ai for gi <C 1. Since C%{E) fluctuates from sample to 
sample, let us make more precise what is the meaning of 
Eq.(2). For the one particle case, the curvature Ci(E) 
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is characterized, in the metallic regime, by a very broad 
distribution (generalized Lorentzian) 113. Its mean value 
over the statistical ensemble is zero at the center of the 
band (E = 0) so that g\ is typically given by the width of 
the distribution, defined for instance by the well behaved 
quantity < |Ci(£')| >. C2(E) should be understood in 
a similar way i.e. as a typical value characterizing the 
width of a distribution, ami not as a mean value. Pre- 
liminary numerical resultaill display a similar behavior 
for the TIP spectrum for a large enough U. 

In the metallic regime, the typical values of C 2 [a) {E) 
and I are given respectively by gi-^e and gi^f so that 
the total curvature is a decreasing function of the inter- 
action through the term oc g 2 gi as expected in a good 
metal. In the insulating regime (gi < 1), C 2 (0) {E) is 
zero up to exponentially small terms while the typical 
value of I is proportional to — Ai so that C 2 increases 
with g 2 in agreement with the derealization effect de- 
scribed above a. 

Although the physical meaning of C 2 is clear, this 
quantity, unlike the case of non interacting electrons in 
a random potential, can hardly be considered to be rele- 
vant for a scaling theory of the interacting system. From 
the random matrix theory point of view, the spectrum of 
the interacting system displays a very unusual behaviour. 
Usually, the curvatures generated by varying some sys- 
tem parameter also set the energy scale below which the 
Wigner Dyson rigidity of the spectrum is observed. This 
is a central point in the scaling theory of localization. 
The fact that g 2 increases with the interaction in the 
metallic regime discards C 2 as a candidate for scaling. 
In other words, we are looking for another "curvature" 
which would be proportional to g 2 (i.e zero for the non 
interacting system). 

Without interaction (U = 0), the TIP-Hamiltonian is 

separable into two identical one particle Hamiltonians, 

2 

i.e. H = hi + h 2 with hi = ^jj + V(xi) and where V(x) 
is the random potential. In other words, the two particles 
ring can be thought of as two independent one particle 
rings. For this argument to be valid, we assume either 
discernable particles, or indiscernable particles with ad- 
ditionnal quantum numbers. This leads us to associate 
distinct Aharonov-Bohm fluxes 0i and <j) 2 to each ring. 
The separabilty of the TIP-Hamiltonian is broken by the 
interaction U(xi—x 2 ). A variation of 0i will then induce 
a current in the second ring characterized by a mutual 
inductance. One can extend this point of view to the case 
where a current is driven due to the interactions between 
two systems, each of them characterized by a well denned 
gauge field (two SNS junctions for instance). Similarly, 
this two particle problem in Id can be thought of as a 
single particle problem in 2d, where by closing the 2d 
system on itself as a torus, one can again introduce two 
Aharonov-Bohm fluxes 0i and <p 2 . Without interaction, 
the 2d character of this equivalent one particle model is 
misleading since the Hamiltonian is separable. But with 
interactions, we have a genuine 2d system and a topology 
similar to those considered, in the description of the inte- 
ger quantum Hall effect E3. We thus define the spectral 



two-form: 
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where \A > are TIP-eigenstates of energies E A . Using the 
structure of the TIP-states in the Fock basis, we obtain: 



C( 01,02) = 02^(01, 02 ) 



(4) 



where U(0i,02) is the Berry connexion whose flux 
through juclosed surface in the flux space is a geomet- 
ric phaseEE This relation constitutes the second result 
of this work. The remainder of this letter is devoted to 
establish and discuss them. 

We consider two particles of mass M interacting 
through a short range attractive or repulsive potential 
U(x\ — x 2 ) and submitted to a random potential V(x). 
Defining, 



M 



C 2 (E)^-J2S(E~E A )- 
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< A\P\A' > 
E A - E' 



(5) 



where P — p\ + p 2 is the total momentum of the two 
particles, the curvature C 2 (E) defined by Eq.(l) rewrites: 



C 2 {E) = B[J2S(E-E A ) + C 2 (E)] 



(6) 



We define the one particle states by (la >,e Q ) and the 
Fock states by (|a/3 >, t a p = e a +£p)- To calculate C 2 (E) 
for the interacting and disordered case, we expand the 
TIP-states \A > on the Fock states la/3 > : 
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so that C 2 (E) contains a product of four complex am- 
plitudes C~j. To proceed further, we assume as is usual 
for disordered metals, that only those trajectories which 
do correspond to time reversed amplitudes contribute 
to the sum in C 2 (E). Therefore, we keep only terms 



of the form 

a = 7,/3 = 



a 



a/3 



and 
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i.e those obtained for 



Combinations of 



the type 8ppi < a! \p\ a > 8 aa i < f3'\p 2 \/3 > cancel since 
< a|pi|a >= for zero magnetic flux. Then, C 2 (E) 
rewrites: 
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C2(E) = -^2S(E-E A )—- 



AA' 



E A - E A > 



££|C , Q/3 ' 4 | \C a >p> A I 5pp\ < a'\pi\a > 

a/3 a'fS' 



In the previous expression appears the density of states 
p a p(E) = EfiS(E -E A )\C^\ 2 For the TIP-problem 
that it is well described by the Breit- 



it was shown 
Wigner form: 
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The shift Tq is negligible for weak enough U and the 
width r can be estimated using the Fermi Golden rule: 
r = 2n < H% d > /A| where < H% d > is the variance of 
the non diagonal matrix elements of H in the Fock basis. 
1 / A| is the effective density of Fock states coupled by the 
two-body interaction. For a lattice model with N sites 
it is assumed that H d are independant normal variables 
characterized by a variance of the order of U 2 /N 3 . This 
estimate comes from the assumption that the one par- 
ticle wavefunctions are ergodic, as implied by the O(N) 
invariance in random matrix theory. This corresponds to 
the zero mode contribution oLa diffusion process (higher 
modes have been considered u) . Within the zero-mode 
approximation, A| = A 2 = A^/B, for gi <C 1. Using 
these approximations we obtain, 

r rm nM it? \ sr I < a '\pi\ a > I 2 

C 2 (E) = — 2_^p BW (E - e a p) ■ 
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The first term in C 2 (E) contributes to C 2 - \E) 1 i.e to the 
two-particule curvature in the abscence of interactions. 
The second part is proportional to the interactions so 
that C2{E) rewrites: 



C 2 (E) 



V ' 2 Ai 
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which defines both C 2 ^ an d I- The sign of C 2 (E) may 
fluctuate from sample to sample but not the relative sign 
between the two terms of the rhs of Eq.(lO). The total 
two-particle curvature is: 



C 2 {E) = C^\E)-l^I 



(10) 



Using A2 = Al. and g 2 = we obtain the Eq.(2) for 
C 2 (E). As discussed in the introduction, the mean value 
of the curvature C 2 is zero. We are then interested in 
either the typical value (when it exists) or in the width 
of the distribution. For the non interacting case (r = 0), 
the two-particle energies are given by Ea — e Q /3 and the 

" d 2 E A (4>) 1 



typical value of the curvature I 



that C 2 (0) 



giAi 



l</.=0 



oc giAi so 



In the same way, replacing J2 a by 



A 1 1 j B B dx, the typical width of I defined by its absolute 
value is: 



Ltyp 
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—1(0, 00). 



(11) 



The constant 1(0, 00) is given by the integral 

1 f°° . 1 1 
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(x + y) 2 + lj_ 00 (z + yf + l \z + y\ 

(12) 



where It yp is evaluated at the center of the band E = 
and in the limit B/T — * 00 where the bandwidth is much 
larger than T. In this metallic limit g\ 3> 1, we obtain 
for the typical two particle curvature: 
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Because of interactions, C 2 is reduced by &, small correc- 
tion. This agrees with numerical results obtained for 
N spinless fermions. Since C 2 (E) can be interpreted as 
a measure of the transport along the system, we recover 
that repulsive interactions decrease the conductivity as 
expected when g\ 3> 1. It could have been anticipated 
from the behavior of the TIP-energy levels Ea(4>)- In the 
absence of interaction (U = T = 0), the TIP-spectrum 
for gi 3> 1 has many level crossings. This results from the 
superposition of two independent spectra (although each 
of them follows the Wigner Dyson statistics). A finite 
interaction U removes these level crossings and there- 
fore reduces their curvature. In contrast, when g\ <C 1, 
level crossings are already suppressed by one particle lo- 
calization and a finite interaction may enhanqe-,the flux 
dependence of Ea (4>) , as numerically observedt-1. 

The relation (11) has been obtained by assuming a uni- 
form broadening of all the Fock states. This is no longer 
correct when g\ < 1, for scales L\ < L < L 2 . There, 
among the N Fock states, only a fraction iV;/(~ LL\ in 
lei) is broadened by the interaction. The remaining part 
Nq m N ~ Njj has a negligible broadening for a short 
range interaction since they do correspond to one parti- 
cle states localized far away from each other. Then, in 
order to calculate C 2 (E) for L\ < L < L 2 , we split the 
sum in Eq.(5) into two sums corresponding respectively 
to Ao TIP-states close to a single Fock state and Njj 
TIP-states broadened over 172 Fock states. The first sum 
corresponds to states \Aq > which are localized within 
L\ <C L. In the Hilbert space of these states, it is always 
possible to build a well defined position operator 113 X 
even for a ring geometry. Therefore, for these states the 



f-sum rule is fulfilled i.e. 



d^E A (<f>) 
37P 
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= and therefore 



C 2 ( \E) is zero up to terms proportional to exp — (LjL\) 
. Then, in the expression of C 2 (E) remains only the sec- 
ond class of states for which the previous argument does 
not hold since they are coupled by the interaction U and 
extended over a scale L 2 > L. These states will still con- 
tribute to C 2 (E) through the second term in the rhs of 

Eq. (10). Since g\ « (f-sum rule), the analog of Eq. (6) 

1 11 1 2 
for one particle yields for M/2\ < a \pi \a > a typical 

value of order Ai. Using this in Eq. (11) and making 
the replacements A 2 -> A| = B/N v and V ->■ T(A|) 
given by the Golden rule with a density 1/A|, we obtain 
for the Njj interaction-assisted states a typical curvature 
C 2 ~ r ^ 2 e ^ when L\ < L < L 2 , i.e. a mobility in- 
creased by the interaction. For L > L 2 this contribution 
decays itself as exp— (L/L 2 ). This is nothing but the 
result of a decimation of the spectrum where all the Ao 
TIP-states not coupled by the interaction do disappear. 
In this regime, the typical curvature C 2 coincides with 
the TIP-conductance g 2 (L) defined by Imry. 
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Since C2 is proportional to gi only in the localized 
regime i.e. when the one particle mobility is suppressed 
by localization, it cannot be considered as an appropri- 
ate scaling parameter for the TIP-localization transition. 
For L <C L\ ( metallic regime), the number 172 of consec- 
utive TIP-levels exhibiting the universal Wigner-Dyson 
rigidity increases with U , but remains much smaller than 
C2 which is essentially dominated by the one particle ki- 
netic part <?i In order to get rid of this one parti- 
cle contribution, we are led to introduce the response to 
two independent fluxes 4>\ and <p2- The simplest gen- 
eralization of a curvature which does not contain the 
one particle kinetic part is the crossed product (two- 
form) defined by Eq.(3). For U — 0, it is straightfor- 
ward to show that C = 0, since \A >= \a[3 > and 

< Wl\ a >< P\-§k > is reaL This is not truc 
anymore for the interacting case. There, using Eq.(7) we 
have: 

C = ImJ2 E C a p A *C je A < ^L| 7 X P\ A > . 



(14) 



where C A p 



-iT((Ea ~ £0/3) — *r) 1 (this results from 
the Lippman-Schwinger equation). The approximation 
previously considered for the C A p i.e. keeping only the 
combinations a — 7 and /3 = e, for the calculation of C2, 
gives a real term which does not contribute. However, by 
taking a = e and /3 = 7, and approximating T((E A — 



e Q /3)) 2 + T 2 ) -1 by A 2 1 , we obtain: 



C=-JmV 

A 2 ^ 

a/3 



dH 



0><l3\ma> 



(15) 



where we recognize the expression of the two-form con- 
nexion V(<j>i, ^2) (see Eq.(4)) whose integral over a closed 
surface in the flux space is a geometric phase which is in- 
dependent of the interaction. 

Similar curvatures of topological origin were consid- 
ered in various contexts. To describe the Mott transi- 
tion in the one dimensionnal repulsive U > Hubbard 
model, Shastry and SutherlandliU introduced two "fluxes" 
coupled respectively to the charge and spin degrees of 
freedom, for a finite density of fermions (no disorder). 
The two flux curvature measures in this case the dif- 
ference between the charge and spin susceptibilities, fa, 
the context of semiclassical physics, Robbins and Berryta 
did consider the semiclassical approximation of a spectral 
two- form similar to those given by Eq.(3) for situations 
where the classical limit corresponds either to integrable 
or chaotic systems. In the integrable case, the semiclas.- 
sical behaviour is equivalent to the Hannay two-forrrO. 
For non integrable systems, physical realizations of such 
a two-form is a debated question. The curvature C we 
introduced for the TIP-problem might be such an exam- 
ple. 
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